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ABSTRACT 

The enhancement factor of the resonant thermonuclear reaction rates is cal- 
culated for the extremely dense stellar plasmas in the liquid phase. In order to 
calculate the enhancement factor we use the screening potential which is deduced 
from the numerical experiment of the classical one-component plasma. It is found 
that the enhancement is tremendous for white dwarf densities if the 12 C + 12 C 
fusion cross sections show resonant behavior in the astrophysical energy range. 
We summarize our numerical results by accurate analytic fitting formulae. 

Subject headings: dense matter — nuclear reactions, nucleosynthesis, abundances 
— plasmas 



1. INTRODUCTION 

In a recent important paper Cussons, Langanke, & Liolios (2002) have pointed out the 
potential resonant screening effects on stellar 12 C + 12 C reaction rates. The 12 C + 12 C fusion 
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cross sections show noticeable resonant structures down to the lowest energies measured so 
far in the laboratory E ~ 2.4 MeV (Kettner, Lorenz-Wirzba, & Rolfs 1980). If the resonant 
structure continues to even lower energies and the astrophysical reaction rate is due to the 
contributions of narrow resonances, one then has to consider that the entrance channel width 
of these resonances will be modified in the plasma. 

Cussons, Langanke, & Liolios (2002) have specifically pointed out the possible impor- 
tance of the plasma effects on the resonant 12 C + 12 C reactions for a carbon white dwarf 
environment with T = 5 x 10 7 K and p = 2 x 10 9 g cm -3 . They have considered a resonance 
energy interval 0.4—2 MeV. They have specifically discussed a rather extreme case of the low 
resonance energy E r = 400 keV and have estimated the overall enhancement of the resonant 
12 C + 12 C reaction rates due to the plasma effects for this case. 

Cussons, Langanke, & Liolios (2002) adopted the method of Salpeter & Van Horn 
(1969) which is based on the lattice model of the dense plasma to calculate the resonant 
screening effects. One of the present authors (N. I.) and his collaborators have calculated 
the enhancement of non-resonant thermonuclear reaction rates in extremely dense stellar 
plasmas (Itoh, Kuwashima, & Munakata 1990). This work is a natural extension of the 
works of Itoh, Totsuji, & Ichimaru (1977) and Itoh et al. (1979), and improves upon the 
accuracy of the results of Salpeter & Van Horn (1969). Itoh, Kuwashima, & Munakata (1990) 
have summarized their numerical results by an accurate analytical fitting formula which will 
be readily implemented in the stellar evolution computations. 

The aim of the present paper is to extend the work of Itoh, Kuwashima, & Munakata 
(1990) to the case of resonant reactions. The present paper is organized as follows. Physical 
conditions relevant to the present calculation are made explicit in § 2. Calculation of the 
enhancement factor of the resonant thermonuclear reaction rates is summarized in § 3. The 
results are presented in § 4. Extension to the case of ionic mixtures is made in § 5. Concluding 
remarks are given in § 6. 



First we consider thermonuclear reactions which take place in the plasma in the thermo- 
dynamic equilibrium at temperature T composed of one kind of atomic nuclei and electrons 
with number densities and n e respectively; Ze and M denote the electric charge and the 
mass of such an ion. The conventional parameters which characterize such a plasma are 
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PHYSICAL CONDITIONS 



r 




(2.1) 
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27vr 2 \ M(Ze) 



1/3 



4 J h?k B T ' (2 ' 2) 

where a is the ion-sphere radius a = (47rnj/3)~ 1 / 3 , A is the mass number of the nucleus, 
T 8 is the temperature in units of 10 s K, and p 6 is the mass density in units of 10 6 g cm -3 . 
In this paper we restrict ourselves to the case that electrons are strongly degenerate. This 
condition is expressed as 



T < T F = 5.930 x 10 9 



{l + 1.018(Z/A) 2 / 3 pf } 



1/2 



K, 



(2.3) 



where Tp is the electron Fermi temperature. Furthermore, we consider the case that the ions 
can be treated approximately as classical particles. The corresponding condition is written 

as 

nA 3 < 1 , (2.4) 

where A = (2nh 2 / MltsT) 1 / 2 is the thermal de Broglie wave length of the ions. The condition 
(2.4) is rewritten as 

T 8 > 2.173pg /3 A4 5/3 , (2.5) 

where p$ is the mass density in units of 10 9 g cm -3 . In this paper we impose a further 
condition that the ions are in the liquid state (Slattery, Doolen, & Dewitt 1982): 

T<178. (2.6) 

The parameter 3r/r corresponds to the ratio of the classical turning point radius at the 
Gamow peak and the mean interionic distance in the case of the pure Coulomb potential 
(Itoh, Totsuji, & Ichimaru 1977; Alastuey & Jancovici 1978). The theories of Itoh, Totsuji, & 
Ichimaru (1977), Itoh et al. (1979) and also Alastuey & Jancovici (1978) are valid under the 
condition 3r/r < 1.6. Itoh, Kuwashima, & Munakata (1990) have extended the calculation 
of the enhancement of the non-resonant thermonuclear reaction rates to the case 3r/r < 5.4. 
In Figure 1 we show the density-temperature diagram of the pure 12 C plasma which illustrates 
the physical conditions described in this section. 



3. ENHANCEMENT OF THE RESONANT THERMONUCLEAR 

REACTION RATES 



In this paper we use the following screening potential for the classical one-component 
plasma 



2^2 



H(r) 



Ck B T - 



Z 2 e 



< r < r 
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r < r < 1.60a 



(3.3) 



(3.2) 



This screening potential is derived from the equilibrium pair correlation function of the 
classical one-component plasma. In adopting this screening potential, our standpoint is the 
same as that of Alastuey & Jancovici (1978). They have argued that the pair correlation 
function should be taken as the static one. The point is that the transmission coefficient of 
the potential barrier is exceedingly small, which makes nuclear reactions very rare events. In 
a loose classical analogy, one might say that, in most collisions, the colliding particles tunnel 
through only a certain distance and are then reflected back. Therefore, as soon as r is larger 
than a few nuclear diameters, equilibrium is achieved and the probability of finding two 
nuclei at a distance r from one another is given by the equilibrium pair correlation function. 
Thus one can use the averaged potential in describing the tunnelling. More recently DeWitt 
(1994), Rosenfeld (1994), and Isern & Hernanz (1994) have studied this problem. They have 
essentially confirmed the correctness of the method of Alastuey & Jancovici (1978) on which 
our work is based. A similar work in the case of the solar fusion reactions has been recently 
carried out by Bahcall et al. (2002) essentially confirming the soundness of the method of 
the average potential. 

We also remark in relaion to the above-stated point that the nuclear reactions as a whole 
are taking place on the macroscopic time scale, whereas the screening potential is kept in 
equilibrium on the microscopic time scale. 

We here remark on the validity of the classical one-component plasma (OCP) model. In 
the interior of dense stars the electron Fermi energy becomes much larger than the Coulomb 
interaction energy between the electron and the ion. Therefore, the electron liquid becomes 
an almost uniform liquid. Owing to this fact, the interior of dense stars can be satisfactorily 
described by the OCP model which consists of point ions imbedded in the rigid background 
of electrons. Of course, at the same density the OCP model becomes better for the smaller 
ionic charge Z, as the ratio of the electron-ion interaction energy to the electron Fermi energy 
becomes smaller for the smaller Z. 

The expression for C is taken from Alastuey and Jancovici (1978). The two segments 
of the screening potential are matched at r = r$ so that the screening potential and its 
derivative be a continuous function with respect to the distance r. This procedure produces 
solutions for r and b for the range of T-values 4 < T < 90. Outside this range we use the 
value of b which makes (3.1) and (3.2) continuous at r—1. 171875a. In this case the first 
derivatives of (3.1) and (3.2) are slightly discontinuous at this point. The linearly decreasing 
part of the screening potential is identical to that employed by Itoh, Totsuji, & Ichimaru 
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(1977) and also by Itoh et al. (1979). The screening potential (3.1) and (3.2) fits the results 
of the numerical experiments excellently. (See Figure 1 of Itoh et al. (1979) for the accuracy 
of this screening potential in reproducing the results of the Monte Carlo computations.) Note 
that this screening potential exactly cancels the Coulomb potential Z 2 e 2 /r at r = 1.60a. We 
further assume that the potential of mean force vanishes for r > 1.60a. Given the explicit 
form of the screening potential, we are now in a position to calculate the enhancement of 
the resonant thermonuclear reaction rates. 

A single resonance in the cross section of a nuclear reaction + 1 — > 2 + 3 can be 
represented most simply as a function of energy in terms of the classical Breit-Wigner formula 
(Fowler, Caughlan, & Zimmerman 1967) 

nh 2 uj r T 1 T 2 

° = (3 - 4) 

^ (£ _ Er )2 + itot 

where [i is the reduced mass, E is the center-of-mass energy, E r is the resonance energy, u r 
is the statistical weight factor, Ti is the partial width for the decay of the resonant state by 
reemission of + 1, T 2 is the partial width for emission of 2 + 3, T tot = Ti + T 2 + • • • is the 
sum over all partial widths. The partial width T\ is proportional to the barrier penetration 
factor P(E) for the screened Coulomb potential. 

Tx oc P{E) = exp j-^^^ r 'V(r)-£] 1/2 drj , (3.5) 
Z 2 e 2 

V(r) = -!L-H(r), (3.6) 
where r tp is the classical turning point radius which satisfies the condition 

V(r tp )-E = 0. (3.7) 

We consider the case that the resonance is sharp; that is, the full width at resonance, 
T r , is considerably smaller than the effective spread in energy of the interacting particles. 
We further consider the case Ti <C T 2 , T tot ~ T 2 . Cussons, Langanke, & Liolios (2002) have 
pointed out that for 12 C + 12 C resonances far below the height of the Coulomb barrier, the 
entrance channel width Ti is much smaller than the total resonance width. The latter (which 
is of order ~ 100 keV for the observed resonances above 2.4 MeV) is also noticeably smaller 
than the resonance energy. In this case we have (Fowler, Caughlan, & Zimmerman 1967) 

. . f 27th 2 \ 3/2 (Uj^r ( E r \ 

M = — — ^^ exp -— ^ , (3.8) 



fik B T J h \ ksT 

(c^7) r = ov7 r = ( ^ lF2 J ~ (ujF^r. (3.9) 

V 1 tot J r . 
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Therefore the partial width r\ in equation (3.5) is to be evaluated at the resonance energy 
E = E r . Here we notice that the resonance energy is shifted by the plasma effects. We take 
E r to be the shifted resonance energy. The shifted resonance energy E r is related to the 
resonance energy in the vacuum E® by the relationship 



E,,. 



E° r - Ck B T, 



(3.10) 



where the expression for C is given by equation (3.3). 
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The barrier penetration factor Pq(E) for the pure Coulomb potential Z 2 e 2 /r of the 
identical nuclei is known to be 



(3.11) 
(3.12) 

Therefore, the enhancement factor a of the resonant thermonuclear reaction rates which 
arises because of the plasma screening effects is 

1 1/2 

71 1 



a = exp 



h 2 



exp(C) , 



j(r» 

h(x) 



rx tp 

Jo 



X 



— h(x) — e 



1/2 



dx , 



a 



H(r) 
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,0 



Z 2 e 2 
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aE° 

Z 2 e 2 ' 



a 



aE r 



r 



(3.13) 

(3.14) 

(3.15) 
(3.16) 

(3.17) 



Here we notice that our method is valid for e r = e° — ^ > 0. 



Since we have (Itoh et al. 2002) 

a = 0.7346 x 10" 
we can rewrite equation (3.13) as 



10 



p 6 \-V3 



P6\ 

a) 



cm , 



(3.18) 



a 



= exp \ -1.004 x 10 1 ZA 2/3 p~ 1/6 



+ </(r> r 



exp(C) 
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exp 



-1.004 x \&ZA*I 3 P - 1/6 K(T, e r )] exp(C) . 



(3.19) 



We also have a useful relationship 

Z 2 e 2 



1.960 x 10 -3 Z 2 



P6 \l/3 



P6\ 
AJ 



MeV. 



(3.20) 



This gives Z 2 e 2 /a = 0.308 MeV and 1.004xl0 1 ZA 2 / 3 pe 1/6 = 99.85 for Z = 6, A = 12, 



p = 10 9 g cm 3 



4. RESULTS 

We have carried out the numerical integration of J(T,e) in equation (3.14) for various 
values of T and e. In Figure 2 we show the function J(T, e) as a function of e for various 
values of T. In Figure 3 we show the function K(T, e r ) as a function of e r for various values 

of r. 

In order to facilitate the application of the numerical results obtained in the present 
paper we will present an accurate analytic fitting formula for K(T, e r ). We have carried out 
the numerical calculations of K(T, e r ) for 1 < T < 200, < e r < 10. We express the analytic 
fitting formula by 

10 

log 10 K(F, e r ) = J2 a ^9 l u\ (4-1) 

i,j=0 

9 = TJ505 (1 ° gl ° r " L15 ° 5) ' (42) 
u = ^ (e r - 5.0) . (4.3) 

The coefficients are presented in Table 1. The accuracy of the fitting is generally better 
than 0.1%. 



5. IONIC MIXTURES 

Itoh et al. (1979) analyzed the results of the Monte Carlo computations for the screening 
potentials of the ionic mixtures of various charge ratios and concentration ratios carried out 
by them and also by Hansen, Torrie, & Vieillefosse (1977). They have established that the 
screening potential at intermediate distances Hij(r) for a mixture of two kinds of ions with 
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charges and number densities (Zie, n\) and (Z 2 e, n 2 ) given below fits the Monte Carlo results 
excellently within the inherent Monte Carlo noise: 



r 



= 1.25 - 0.390625 



(a» + CLj)/2 



ZiZje 2 
(l/2){ ai + aj )k B T 

3Zi 

_47r(Zini + Z 2 n 2 ) 
3Z 2 

Air(Z 1 n 1 + Z 2 n 2 ) 
We also define the parameter 

'27tt 2 \ 2//, ; ZfZ- ( 



Oi = 



a 2 = 



= 1,2), 
(i,j = l,2), 



1/3 



1/3 



(5.1) 
(5.2) 

(5.3) 

(5.4) 



(z,j = l,2), 



(5.5) 



where //y is the reduced mass for the two ions with charges Zj and Zj. Then the enhancement 
factor for the resonant thermonuclear rates of the two nuclei Zi and Zj is given by 



a = 



= exp < 


-2 


-K + ai) 


/I 2 






2nijZiZjC 


= exp < 


f 

-2 


' 1 r ^' 

-{Oi + aj) 


h 2 






2fiijZiZje 



1/2 



1 



J(Xij, 



ij, t r , 



> exp(Ci 



1/2 



K(T lv e r ) 



exp(C ii ) , 



e,. = - 



1 (oj + Oj)i? r 1 (oj + ctj)-E° Cjj 

2 ZiZj-e 2 ~ 



2 ZiZj-e 2 

1/4 



r ' 



Cij = 1.0531 r i7 - + 2.2931 r*f - 0.5551 In r i7 - - 2.35 . 



y 



(5.6) 

(5.7) 
(5.8) 



The analytic fitting formula for this case if(ry,e r ) has the same form as equations (4.1), 
(4.2), (4.3). 



6. CONCLUDING REMARKS 



We have presented a calculation of the enhancement of the resonant thermonuclear 
reaction rates for extremely dense stellar plasmas. The calculation has been carried out by 
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adopting the screening potential derived from the Monte Carlo computations of the classical 
one-component plasma. We have summarized our numerical results by an accurate analytic 
fitting formula to facilitate applications. The present results will be useful if the 12 C + 12 C 
fusion reaction contains narrow resonances in the astrophysical energy range. 

We wish to thank K. Langanke for making the preprint available to us prior to its 
publication. We also thank Y. Oyanagi for allowing us to use the least-squares fitting program 
SALS. We are most grateful to our referee for many valuable comments on the original 
manuscript which helped us tremendously in revising the manuscript. This work is financially 
supported in part by Grants-in-Aid of the Japanese Ministry of Education, Culture, Sports, 
Science, and Technology under contracts 13640245, 13740129. 
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Table 1 


: Coefficients 


(Xj -) 












7=0 


7=1 

J 


7=2 


7=3 


7=4 


7=5 


i=0 


-3.65927E+0 


-1.26734E+0 


6.07584E-1 


6.30602E-2 


-7.78163E-1 


-1.79849E+0 


i=l 


1.55574E-3 


-4.55740E-3 


2.58322E-2 


7.00779E-2 


-2.64033E-1 


-3.01686E-1 


i=2 


6.20281E-3 


-1.79841E-2 


1.04425E-1 


2.79543E-1 


-1.07405E+0 


-1.21880E+0 


i=3 


-7.53523E-4 


3.65732E-2 


-2.62102E-1 


-7.31969E-1 


2.74892E+0 


3.19750E+0 


i=A 


1 18877E— 2 


1 11Q7QE— 2 


— 1 051 ORE— 1 


—4 62001 E— 1 


1 32482E+0 


2 26167E+0 


i=5 


-4.10209E-2 


-9.73313E-2 


8.96704E-1 


2.55956E+0 


-9.59895E+0 


-1.13011E+1 


i=6 


1.89159E-2 


5.88188E-2 


-5.47015E-1 


-1.22948E+0 


5.38904E+0 


4.80222E+0 


i=7 


5.10796E-2 


1.10523E-1 


-1.09263E+0 


-2.99190E+0 


1.16478E+1 


1.32623E+1 


i=8 


-4.26759E-2 


-1.08244E-1 


1.02972E+0 


2.53992E+0 


-1.04973E+1 


-1.05685E+1 


i=9 


-1.97340E-2 


-4.38567E-2 


4.51951E-1 


1.19137E+0 


-4.79595E+0 


-5.30281E+0 


z=10 


2.00491E-2 


5.10687E-2 


-4.93422E-1 


-1.22559E+0 


5.06476E+0 


5.18690E+0 


















j=6 


M 


.7=8 


.7=9 


.7=10 




i=0 


3.53738E+0 


2.90085E+0 


-4.98720E+0 


-1.80622E+0 


2.62678E+0 




i=l 


8.84245E-1 


4.72177E-1 


-1.18370E+0 


-2.24453E-1 


5.27141E-1 




i=2 


3.62099E+0 


1.92629E+0 


-4.87413E+0 


-9.25366E-1 


2.18368E+0 




i=3 


-9.25572E+0 


-5.12232E+0 


1.25061E+1 


2.53278E+0 


-5.66738E+0 




i=4 


-4.95541E+0 


-4.03418E+0 


7.40638E+0 


2.26404E+0 


-3.72743E+0 




i=5 


3.24488E+1 


1.84055E+1 


-4.41392E+1 


-9.35145E+0 


2.02577E+1 




i=6 


-1.70172E+1 


-6.83208E+0 


2.14258E+1 


2.91262E+0 


-8.99871E+0 




i=7 


-3.92927E+1 


-2.16796E+1 


5.33707E+1 


1.10456E+1 


-2.44789E+1 




i=8 


3.41813E+1 


1.61786E+1 


-4.46244E+1 


-7.64260E+0 


1.95906E+1 




i=9 


1.61526E+1 


8.69953E+0 


-2.19164E+1 


-4.44246E+0 


1.00465E+1 




?;=io 


-1.66155E+1 


-8.08293E+0 


2.18847E+1 


3.90128E+0 


-9.71024E+0 
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Fig. 1. — Physical conditions for the pure 12 C plasma. The present calculation is valid in the 
region bounded by the electron Fermi temperature Tp and the line r^A 3 = 1. Furthermore, 
the ions should be in the liquid state: T < 178. The lines of r=const., T=const., and 
3r/r=const. are shown. 
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Fig. 2. — J(r, e) as a function of e for various values of T. 
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SIS db function of e r for various values of T. 



